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Subsurface porous formations provide large capacities for underground hydrogen storage (UHS). Successful uti- 
lization of these porous reservoirs for UHS depends on accurate quantification of the hydrogen transport char- 
acteristics at continuum (macro) scale, specially in contact with other reservoir fluids. Relative-permeability and 
capillary-pressure curves are among the macro-scale transport characteristics which play crucial roles in quan- 
tification of the storage capacity and efficiency. For a given rock sample, these functions can be determined if 
pore-scale (micro-scale) surface properties, specially contact angles, are known. For hydrogen/brine/rock system, 
these properties are yet to a large extent unknown. In this study, we characterize the contact angles of hydrogen 
in contact with brine and Bentheimer and Berea sandstones at various pressure, temperature, and brine salinity 
using captive-bubble method. The experiments are conducted close to the in-situ conditions, which resulted in 
water-wet intrinsic contact angles, about 25 to 45 degrees. Moreover, no meaningful correlation was found with 
changing tested parameters. We monitor the bubbles over time and report the average contact angles with their 
minimum and maximum variations. Given rock pore structures, using the contact angles reported in this study, 
one can define relative-permeability and capillary-pressure functions for reservoir-scale simulations and storage 


optimization. 


1. Introduction 


Successful transition towards low-carbon energy systems depends 
not only on harvesting more renewable resources but also on advance- 
ments of large-scale (TWh) storage technologies. Renewable energy can 
be stored in TWh scales, if it is converted into green molecules such as 
hydrogen. The green hydrogen can then be stored in underground ge- 
ological formations, e.g., in depleted hydrocarbon reservoirs and saline 
aquifers (Gabrielli et al., 2020; Stone et al., 2009; Rudolph, 2019). 

Several research studies and a few pilot tests related to underground 
hydrogen storage (UHS) in porous reservoirs have been recently initi- 
ated (Tarkowski, 2019; Panfilov et al., 2012; Panfilov, 2016; Nemati 
et al., 2020; Ebrahimiyekta, 2017). UHS in porous formations still re- 
mains a challenge, due to lack of characterization data needed as in- 
put parameters to perform reservoir simulation and robust storage op- 
timization. Among these input parameters, hydrogen surface properties 
in contact with reservoir fluids, especially brine, are crucially important 
(Heinemann et al., 2021). 

To date, there exists only two experiments for characterizing hydro- 
gen contact angle in subsurface systems. First, a coreflooding test was 
performed in which hydrogen was injected in brine-saturated Vosges 


* Corresponding author. 


sandstone rock at two different pressure (P) and temperature (T) val- 
ues of (50 bar, 20 °C) and (100 bar, 45 °C) (Yekta et al., 2018). It 
resulted in receding contact angles of 21.56° and 34.9°, respectively, 
for the first and second (P,T) values (Yekta et al., 2018). More recently, 
another study was performed in which receding and advancing con- 
tact angles were measured using tilted plate experimental technique 
(Iglauer et al., 2021). That study was performed at pressure range of 
0.1-25 MPa and temperature range of 296-343 K. Both pure and aged 
quartz samples were used with stearic acid in contact with brine (10 wt% 
NaCl). The study reported that the increase of pressure or temperature 
resulted in increasing contact angles from 0° to maximum of around 
50° for pure quartz. However, when the quartz samples were aged for 
several months with stearic acid, intermediate wetting conditions were 
observed (Iglauer et al., 2021). 

Despite its crucial impact in successful development of UHS technol- 
ogy, to date, there exists no hydrogen-specific dataset in general, nor 
contact angle measurements across scales, to allow for reliable site selec- 
tion, development, and storage optimization (Ebigbo et al., 2013; Hage- 
mann et al., 2016; Sainz-Garcia et al., 2017; Tarkowski, 2019; Lubon and 
Tarkowski, 2020; Flesch et al., 2018; Ozarslan, 2012; Panfilov et al., 
2012). Specially once the contact angles are known, for a given rock 
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type, one can perform pore-scale modelling to find upscaled relative- 
permeability and capillary-pressure curves (Hashemi et al., 2021; Klise 
et al., 2016; Liu et al., 2018; Berg et al., 2016). These functions will be 
used as input parameters for reservoir scale simulation studies (Blunt, 
2017; Kunz et al., 2018; Rücker et al., 2019). Note that contact angles in 
cyclic storage transport is hysteretic (Morrow, 1975; Yang et al., 1999). 

To resolve this knowledge gap, in this work, we perform contact an- 
gle measurements for hydrogen/brine/sandstone rock using a captive- 
bubble cell device. Since, there is no external viscous forces to displace 
fluid and gas phases, our study allows for nearly static (intrinsic) con- 
tact angle measurements for hydrogen when it comes in contact with 
the saturated porous reservoir rock. We first benchmark our measure- 
ments for nitrogen gas with the published literature, and then introduce 
hydrogen gas in the system. Table 8 and Table 9, in the appendix, pro- 
vide a summary of different contact angle measurement methods in the 
systems of gas/brine/solid surface. The study is performed under differ- 
ent pressure, temperature and salinity of the brine; so to represent a fair 
assessment of the in situ conditions. Our findings shed new lights on the 
characteristics of the hydrogen surface characteristics when it is stored 
in the subsurface reservoirs. 

The structure of the paper is as follows. First, the experimental setup 
and procedure to measure the contact angle will be described in detail. 
Then the image processing methodology to measure the contact angles 
from the captured images is presented. In its following section, results 
and their relevance for UHS will be discussed. Finally, concluding re- 
marks are presented. 


2. Methods and materials 


In this study, the Captive-Bubble method was utilized to measure 
intrinsic contact angles using a gas bubble at a solid-liquid interface. 
This method is advantageous over the alternative Sessile Drop method, 
because spreading and diffusion of the brine into porous hydrophilic 
substrates in the latter method poses experimental challenges, making 
the data less reliable (Prydatko et al., 2018). 


2.1. Materials 


Hydrogen with a purity of 99.99 mol%, produced by Lindegas, 
was used. The rock slabs were taken from a sawed homogeneous Ben- 
theimer and Berea sandstone blocks. The properties of Bentheimer and 
Berea have been studied in literature (Peksa et al., 2015; Farajzadeh 
et al., 2017; Kapetas et al., 2015). The samples were mainly com- 
posed of quartz (95%)which was evenly distributed throughout the rock 
matrix (see appendix Figs. 10 and 11). The average porosity of the 
Bentheimer and Berea sandstones were about 20%, and permeability 
around 2 to 3, and 0.1 Darcy, respectively. Each slab had dimensions of 
30 x 6 x 12 mm. In addition to brine containing NaCl, a synthetic sea- 
water (Laskaris, 2015) with the composition in Table 1 was also used in 
the experiments. 
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Table 1 

Synthetic seawater composition. 
Salt Quantity (ppm) 
NaCl 24,500 
KCI 670 
MgCl2.6H20 10,150 


CaCl2.2H20 1450 


2.2. Microscopic image analysis of rock samples 


To quantify the surface roughness of rock slabs, 2D and 3D mi- 
croscopic images were taken using a LEICA 3D stereo explorer Fig. 1. 
The surface profiles were characterized based on the internationally- 
recognized standard of EN ISO 287, where the so-called P, factor defines 
the surface roughness (Kaveh et al., 2014). P, is the arithmetic mean of 
the absolute ordinate height values Z(x) within a sampling length (1,). 
The average roughnesses of the Bentheimer and Berea slabs have been 
measured 0.030 and 0.025 mm, respectively. 


2.3. Experimental setup 


The captive-bubble setup, modified after Kaveh et. al. (2014), was 
used to measure contact angles in the system of hydrogen/brine/rock 
at high pressures and temperatures. The schematic of the experimental 
setup is given in Fig. 2. The setup comprises of a single steel cell, which 
holds the rock sample and brine. The injection of the brine and other 
fluids takes place from the bottom and from different inlets. Extraction 
of the brine and gas from the cell was done from the top. The pressure in 
the cell was held constant by a back-pressure regulator connected to a 
large nitrogen cylinder. Injection of both brine and hydrogen was done 
with the use of two Vindum pumps. To bring the brine and hydrogen 
to equilibrium, a pressure gauge was installed in between the hydro- 
gen pump and the cell. Hydrogen flows through a line with respective 
inner and outer diameters of 0.25 mm and 1.58 mm to a nozzle from 
which the hydrogen was released into the brine. The bubble that was 
created on the rock surface was photographed with a Canon 90 camera 
(with a maximum resolution of 12.3 MP) attached to an endoscope. The 
pressure and temperature of the setup were monitored continuously and 
recorded in a computer connected to the system. 


2.4, Experimental procedure 


To start the experiment the cell, including the rock sample, was ini- 
tially filled with the brine over night to reach equilibrium. Then, at 
the desired pressure and temperature, a hydrogen bubble was injected 
through a needle at the bottom of the cell. Due to buoyancy, the bub- 
ble arises and sticks underneath of the rock sample. The pictures of 
the released bubbles were taken by using the connected digital cam- 
era through the endoscope at one side of the cell and diffusion of the 


Fig. 1. Microscopic surface images: a) 2D- 
Bentheimer, b) 2D-Berea, and c) 3D-Berea. 
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light source at the other side. To get sharp pictures of the bubbles, the 
resolution of imaging was set to 6.9 MP (3216 x 2136). Due to diffusion 
and dissolution of hydrogen into brine, the size of the bubbles changed 
over time. Consequently, several images were taken after injection of 
one bubble at each pressure and temperature. Therefore, for each ex- 
perimental condition, the minimum and the maximum contact angles 
were reported in addition to the averaged values. 


2.5. Image analysis 


To calculate the contact angles, the captured images from the hydro- 
gen bubbles were analyzed using an in-house MATLAB code. The needle 
diameter inside the cell was used to define the scale of the images. After- 
wards, the image was converted to grey-scale format and was cropped 
to only keep the bubble and the rock surface. To find the boundary of 
the bubble, the image was subsequently binarized. Tracing boundary of 
the bubble and detecting the apex as well as contact points were fol- 
lowed by fitting the best curve based on the Axisymmetric Drop Shape 
Analysis-Profile (ADSA-P) technique (Li et al., 1992). All the main steps 
of the image analysis are shown in Fig. 4. The ADSA-P technique fits 
the best theoretical Laplacian curve on the physical observed bubble in- 
terface (Li et al., 1992). The Young-Laplace capillarity equation for two 
fluid phases is given as 


1 1 
AP =o(—+—), 1 
o( R; Ro? a) 
where o is interfacial tension, R, and R, are the two principle radii of 
the curvature. Because of the axisymmetry of the bubble in Eqn 1, the 


radii are considered equal at the apex (R; = R3 = Ro), i.e., 


20 
AP apex = Ro (2) 
Also, by considering the gravity as the only external force across the 
interface, the pressure difference in Eqn 2 is assumed to be a linear 
function of the hydrostatic pressure (Apgz) with interception of AP) at 
a reference plane, i.e., 


AP = AP) + Apgz. (3) 


As shown in Fig. 3, the origin of the coordinate system was placed at the 
apex point. The x-axis is tangent to the origin and normal to the axis of 
symmetry. Therefore, Eqn 3 can be rewritten as 

sin(@), _ 20 


=—+A 4 
ae) i” P&Z, (4) 


1 
— + 
Oe 
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Fig. 2. Schematic of the captive-bubble cell exper- 
imental setup. 


Vindum 


Fig. 3. Schematic of an axisymmetric drop, modified after Li et al. (1992). 


where R, rotates in the plane of x-z and Ry = x/sin(@) rotates in the 
plane of z-y. Eqn 4 is restated as 


dx 


T cos0 (5) 
z = sin (6) 


dð _ 2 P Apgz 2 sin (7) 
ds Ro o x 

Eqn 5, 6, and 7 are first order differential equations in terms of the arc 
length (s). Finally, these three differential equations (i.e., Eqn 5, 6, and 
7) are integrated using Runge-Kutta numerical approach to find theo- 
retical Laplacian curve. In order to find the corresponding contact an- 
gle, the objective function is defined to minimize the deviation of the 
physically-observed curve from a theoretical curve by adjusting param- 
eters: radius of the curvature at the apex, Rg and interfacial tension, o. 
The detail of this procedure has been extensively explained in the lit- 
erature (Li et al., 1992). Brine and gas density were calculated based 
on the literature thermodynamic formulation for each specific pressure 
and temperature condition (Batzle and Wang, 1992) and reported in the 
appendix, Tables 10-15. 


2.6. Calibration of the setup 


Before starting the main test for the hydrogen gas, the setup was 
calibrated against the literature data for a nitrogen/brine/quartz sys- 
tem. The literature reports contact angle of nitrogen on a smooth alpha- 
quartz crystal surface at pressure of 13 MPa and temperature of 333K 
(Al-Yaseri et al., 2015). Under the same experimental conditions, the 
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Table 2 
Summary of calibration of the setup with the literature data (Al-Yaseri et al., 2015) for nitro- 
gen/brine/sandstone. 

Parameters Literature test (Al-Yaseri et al., 2015) Our test 


Liquid phase 5000 ppm NaCl 


5000 ppm NaCl 


Gas phase N3 N3 

Rock sample Smooth alpha-quartz crystal surface Bentheimer (~95% quartz) 
Pressure (bar) 130 130 

Temperature (°C) 60 60 

Contact angle (°) 40.6 +3.9 40.8 +5.9 


method used in our study (i.e., captive-bubble method), resulted in 
the contact angle of (40.8° +5.9°) for nitrogen/brine/Bentheimer sand- 
stone system, which is in close agreement with the reported data (Al- 
Yaseri et al., 2015), i.e, 40.6° +3.9°. Summary of the validation test is 
provided in Table 2. 


2.7. Summary of the test cases 

Extensive experiments were conducted to examine the effect of pres- 
sure, temperature, salinity, and rock type on the contact angle in the hy- 
drogen/brine/sandstone system. The experimental conditions are sum- 
marized in Fig. 5. 
3. Results and discussion 

In this section the results of the experiments will be explained in 
detail. The impact of each parameter, i.e., pressure, temperature, salin- 
ity, rock type will be separately discussed. Note that, in order to prove 
reproducibility of the results, measurements were repeated up to three 
times (see appendix Table 16-17--18). 
3.1. Effect of bubble size on the contact angle 


As shown in Fig. 6, it was noticed that due to dissolution/diffusion 
of the hydrogen gas into brine, the size of the created bubbles contin- 
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b) Gray-scale image 
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uously decreased during the experiments until they disappeared. Inter- 
estingly, as the size of the bubble decreases, the calculated contact angle 
increases. A similar behavior has been reported by Kaveh et. al. (2014), 
Haeri et al. (2020), and Jung et al., (2012) for CO,/brine/rock system. 
Therefore, to capture the effect of the bubble size on the reported contact 
angles, for each test case, several images were taken from every injected 
bubble at different times, out of which only four images were analyzed 
to calculate the minimum, maximum, and the mean contact angle. The 
mean contact angle is the arithmetic average of the four measured con- 
tact angles. An example is given in Fig. 6a, which results in the contact 
angles shown in Fig. 6b. 


3.2. Test Case 1: Effect of pressure and temperature 

Fig. 7 shows the effect of different pressures and temperatures on the 
contact angles of the hydrogen/water/Bentheimer system in the absence 
and presence of NaCl (5000ppm). No obvious correlation was found, 
as all the data points fall within the accuracy range of the conducted 
experiments. Detailed results are given in Tables 3—4. 
3.3. Test Case 2: Effect of salinity 

To quantify the impact of salinity, brines with three different salin- 


ities, i.e., 0, 5000 and 50,000 ppm NaCl were used at a constant 
temperature of 30 °C and four different pressures in the range of 20 


Fig. 4. The image processing procedure. 


1.5 


e) Fitted bubble profile by using Young-Laplace equation 
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Table 3 
Contact angle values of hydrogen/pure water/Bentheimer. 
Test No. Temp. (°C) Press. (bar) aye (°) Orange (°) Vol.ave (mm?) Vol. ange (mm?) 
1 T~20°C 22.3 20.3 30 [28.7, 32.8] 4.56 [2.36, 7.05] 
23.5 50.2 32.6 [29.2, 39.1] 5.49 [1.34, 9.96] 
23.4 70.7 31.1 [25.9, 37.8] 3.93 [1.45, 9.10] 
23.9 100.5 30 [26.0, 36.9] 4.05 [1.48, 7.36] 
2 T~30°C 31.9 22 33.7 [30.6, 37.1] 3.48 [2.21, 4.66] 
32.5 51.8 30.5 [29.4, 32.9] 3.09 [2.20, 3.66] 
32.8 71.5 33.9 [32.6, 36.5] 3.39 [2.38, 4.48] 
33.2 100.5 31.7 [29.0, 39.0] 5.27 [1.93, 9.49] 
3 T~40°C 39.5 20.3 31.9 [29.0, 35.3] 5.09 [2.61, 8.38] 
39.9 50.2 29.8 [26.3, 35.9] 7.42 (2.21, 12.67] 
40.1 72.8 31.2 [28.9, 36.0] 7.04 (2.27, 12.58] 
40.3 100.3 32 [28.9, 35.2] 3.91 [2.32, 6.02] 
4 T~50°C 49.1 19.8 28.4 [26.1, 29.2] 7.42 [3.96, 10.65] 
49.2 50.6 33.2 [29.4, 39.3] 4.7 [1.68, 8.39] 
49.3 70.2 29.8 [28.6, 31.2] 4.41 [2.66, 6.33] 
49.3 101.2 32.8 [29.9, 38.0] 4.12 [2.14, 6.35] 
Table 4 
Contact angle values of hydrogen/brine (5000 ppm NaCl)/Bentheimer. 
Test No. Temp. (°C) Press. (bar) aye (°) Orange (°) Vol. ve (mm?) Vol. ange (mm?) 
1 T~20°C 21.3 20 33.1 [30.0, 39.2] 3.99 [1.62, 6.14] 
22.1 51.9 29.1 [26.8, 32.8] 4.08 [1.82, 6.53] 
22.3 71.5 33.5 [29.3, 40.5] 3.76 [1.28, 6.50] 
22.9 100.5 33.9 [29.7, 42.7] 4.13 (1.01, 7.37] 
2 T~30°C 38.9 21 29.5 [28.7, 30.5] 4.61 [2.67, 6.55] 
32.2 49.9 34.9 [30.8, 42.2] 3.42 (1.21, 5.77] 
32.7 71.1 36 [32.8, 41.6] 2.8 [1.19, 4.72] 
33.1 98.9 31.9 [31.1, 34.1] 5.59 [2.08, 11.13] 
3 T~40°C 38.9 19.6 32.7 [30.7, 36.0] 4.33 [2.45, 6.32] 
39.5 50.8 34.1 [30.5, 40.7] 3.91 [1.44, 6.48] 
39.9 69.9 34.3 [29.4, 43.0] 4.06 [1.16, 7.18] 
40.1 100.1 37.3 [34.0, 41.3] 2.24 [1.34, 3.22] 
4 T~50°C 47.4 20.7 33.6 [29.2, 40.2] 4.48 [1.51, 7.78] 
48.3 51.3 33.6 [29.9, 41.4] 4.03 [1.40, 6.45] 
49 70.6 34.2 [30.0, 41.5] 4.34 [1.50, 7.91] 
49.2 100.7 33.7 [29.9, 41.6] 5.7 (1.31, 12.66] 
Gas phase Rock phases Brine phases Temperature Pressure 
(°c) (bar) 


Re a S] 20, 30, 40, 50 mss] 20,50, 70, 100 
Sandstone: ee al 20, 30, 40,50 Fe) 20, 50, 70, 100 
Bentheimer 5 0 

ceed a 20, 50, 70, 100 
E a fi ii zi 
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Fig. 5. Experimental conditions for hydrogen tests. 


to 100 bar. The results are shown in Fig. 8. The change of salinity 3.4. Test Case 3: Effect of rock type 

did not result in a meaningful change in the measured contact an- 

gles, indicating that the wetting state of the rock was insensitive to Fig. 9 shows the effect of rock type on the contact angles of the hy- 
salinity in the presence of hydrogen. Detailed results for salinity of drogen/pure water/rock system with Bentheimer and Berea sandstones. 
50,000 ppm NaCl and seawater are also provided in Tables 5 and No obvious correlation was found, as all the data points fall within the 
6. accuracy range of the conducted experiments. Detailed results are given 
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Fig. 6. Effect of bubble size in the system of hydrogen/water/Bentheimer at 23.5°C and 51.2 bar: a) Volume changes over time, b) Corresponding reported range of 


contact angles. 
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Fig. 7. Effect of pressure and temperature in the system of hydrogen/brine/Bentheimer at 20, 30, 40 and 50 °C and 20, 50, 70 and 100 bar, a) pure water and b) 


brine with 5000 ppm NaCl. 


in Tables 3 and 7. To confirm repeatability of the tests, the results of the 
repeated experiments are reported in Tables 16, 17, 18 and 19 in the 
appendix. 


4. Conclusions and outlook 


Wettability of the rock in contact with brine and hydrogen plays 
a crucial role in the displacement processes in UHS. This paper re- 
ports experimental measurements of the contact angle of the hydro- 
gen/brine/sandstone system, relevant for underground hydrogen stor- 
age. This is the first step in understanding and quantifying the impact 
of different parameters in accurate prediction of the fate of stored and 
produced hydrogen. The captive-bubble method was used for measuring 


contact angle. To ensure accurate measurements, the setup was success- 
fully calibrated against the existing literature data with nitrogen gas. 
Then, through several test cases, the intrinsic contact angles were mea- 
sured under various experimental conditions. It was found that, under 
our experimental conditions and within the accuracy of the method 
and setup, the contact angles in the examined systems were not af- 
fected by temperature, pressure and salinity. Under all conditions, the 
presence of hydrogen did not appear to affect the wettability of the 
Bentheimer sandstone. All the results indicated water-wet conditions 
with contact angles in the range of 21.1° to 43°. This general conclu- 
sion and contact angles of less than 50° agree well with the conclu- 
sions of literature (Yekta, Manceau, Gaboreau, Pichavant, Audigane, 
2018). A major discrepancy between our data and the data reported 
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Fig. 8. Effect of salinity in the systems of hydrogen/brine/Bentheimer using pure water, brine (5000 ppm NaCl), brine (50,000 ppm NaCl), and seawater at about 
30 °C for four pressure values of 20, 50, 70, 100 bar. 


Table 5 
Contact angle values of hydrogen/brine (50,000 ppm NaCl)/Bentheimer. 
Test No. Temp. (°C) Press. (bar) aye (°) Orange (°) Vol.ave (mm?) Vol. ange (mm?) 
1 T~30°C 31.3 21.1 333 [30.4, 36.6] 4.02 [2.10, 6.49] 
31.9 51.4 32.8 [30.3, 37.5] 3.54 [1.66, 5.34] 
33 70.6 31.6 [29.1, 36.7] 3.31 [1.53, 5.38] 
33.3 100.7 34.5 [29.8, 42.5] 3.84 [1.36, 5.77] 
Table 6 
Contact angle values of hydrogen/brine (synthetic seawater) /Bentheimer. 
Test No. Temp. (°C) Press. (bar) aye (°) Orange (°) Vol.ave (mm?) Vol. ange (mm?) 
1 T~30°C 31.0 20.5 27.47 [20.8, 35.6] 5.61 [2.62, 8.72] 
31.3 50.2 35.5 [29.1, 43.5] 1.4 [0.42, 2.64] 
31.6 69.6 38.3 [31.7, 44.4] 1.38 [0.61, 2.09] 
31.6 100.7 34.7 [28.8, 41.7] 2.61 [1.11, 4.42] 
Table 7 
Contact angle values of hydrogen/pure water/Berea. 
Test No. Temp. (°C) Press. (bar) aye (°) Orange (°) Vol.ave (mm?) Vol. yange (mm?) 
1 T~20°C 23.6 20.8 30.4 [27.5, 34.9] 4.96 [1.14, 9.19 
23.5 50.6 29 (26.8, 34.5] 5 (0.96, 10.74] 
23.7 70.2 29.1 [25.2, 33.3] 5.22 [1.04, 10.28] 
23.9 100.7 29.6 [23.6, 41.9] 4.17 [0.33, 9.51 
2 T~30°C 32.6 19.4 26.1 [23.3, 29.2] 7.65 [1.36, 15.81] 
32.7 50 23.6 [21.1, 27.9] 7.28 (0.79, 14.17] 
32.8 69.3 31.2 [27.9, 36.1] 3.38 [0.54, 7.29 
33 101.1 31.7 [28.3, 39.3] 3.61 [0.45, 8.59 
3 T~40°C 38.6 21.2 31.1 [27.3, 34.3] 3.01 [0.95, 5.64 
38.6 51 29.5 [25.1, 34.8] 4.34 [0.94, 9.68 
38.6 69.4 29.4 [27.4, 32.5] 3.56 [1.17, 7.14 
38.9 100.7 28.9 [25.8, 31.6] 4.98 (1.52, 10.35] 
4 T~50°C 47.6 20.5 27 [25.6, 30.5] 5.57 [0.91, 10.51] 
47.8 49.4 26.4 [23.2, 31.4] 6.49 (0.93, 14.03] 
48.2 70.6 30.1 [29.4, 31.5] 6.68 [5.13, 8.52 
48.2 99.7 30.5 [27.2, 35.3] 4.55 [0.91, 9.80 ] 


by Iglauer et al. (2021) is the absence of a general trend between the ues due to dominance of gravity for the larger gas bubbles, and possibly 


parameters (pressure, temperature and salinity) and the measured con- different surface composition of the rock. 

tact angles. This can be attributed to the differences in the measurement The presented experimental data and images are all made digi- 
methods, experimental conditions, and sample preparation. It was also tally available open-source at https://gitlab. tudelft.nl/ADMIRE_Public/ 
observed that the bubble size affects the calculated contact angle val- PoreScale_H2 repository. 
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Fig. 9. Effect of rock type in the system of hydrogen/pure water/rock at 20, 30, 40 and 50 °C and 20, 50, 70 and 100 bar, a) Bentheimer and b) Berea. 
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Fig. 10. SEM data of Bentheimer sandstone. 
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Fig. 11. SEM data of Berea sandstone. 


Table 8 
Review of literature experimental measurement of contact angles for systems containing Gas/Brine/Rock. 
Reference Substrate Fluids Gases Method P &T conditions Comments Wettability 
Yang et al. (2008) Carbonate Weyburn Weyburn CO, Sessile-Drop 0.1-31.6 MPa; @ increases with P, intermediate wet 
limestone brine ADSA 300 & 331 K @ decreases with T 27°C (45-100)° 58°C 
limestone (88-135)° 
Espinoza and Carbonate, Calcite water, NaCl, CO, Sessile-Drop up to 20 MPa; @ increases with P, on strongly water-wet 
Santamarina, (2010) Sandstone, Quartz, 200 g/l ADSA 298 K nonwetting surfacs: Calcite: 20° Quartz: 
oil-wet Quartz PTFE PTFE & oil-wet Quartz 20° Oil-wet Quartz: 
e slightly decreases in 85-90° PTFE: 
water-wet Quartz & 100-140° 
Calcite surfaces, 
Ə increases with 
salinity 
Bikkina, (2011) Carbonate, Calcite; water CO, Sessile-Drop 1.3-21 MPa; @ increases with P weakly water-wet 
Sandstone, Quartz 298-323 K slightly then decreases 
slightly and remains 
fairly constant 
Mills et al. (2011) Carbonate, Calcite Brine CO, Nz Captive-Bubble 6 MPa & @ increases with P strongly water-wet 
Sandstone, Quartz 295 K; 13 MPa slightly 
Biotite Orthoclase & 313 K 
Labradorite 
Broseta et al. Carbonate Calcite NaCl, 4g/1 CO, Captive-Bubble 0.5-14 MPa; @ increases with P weakly water-wet 
(2012) 308 K 
Rousse caprock (70% NaCl, 4g/l, CO, Captive-Bubble 1.1-15.5 MPa; e remains constant strongly to weakly 
Calcite) 40g/1, 350¢/1 3.8-413 K with P, ð increases water-wet 
with salinity 
Jung and Sandstone and Silica water, brine CO, Captive- 0.1-25 MPa & @ increases with P in: stronglywater-wet 
Wan (2012) 5M Bubble;Sessile- 318 K;20 MPa 7-10 MPa, after 
Drop;No & 323 K 10 MPa@ remains 
difference in constant @ increases 
results with salinity linearly 
Wang et al. (2013) Carbonate and Calcite water, co, Captive-Bubble 7 MPa & @ increases with P strongly water-wet 
Synthetic 303 K; 20 MPa slightly, @ decreases 
brines & 323 K with salinity 
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Table 9 
Review of literature experimental measurement of contact angles for systems containing Gas/Brine/Rock (continue). 
Reference Substrate Fluids Gases Method P &T conditions Comments Wettability 
Farokhpoor et al. Carbonate Calcite water, NaCl, CO, Captive-Bubble 0.1-40 MPa; @ increases with P on strongly water-wet 
(2013) Sandstone Quartz 0.2 &0.8 M 309 K Mica for the rest not 
feldspar Mica significant changes 
Ə increases with T on 
Quartz 
Kaveh et al. (2014) Sandstone Bentheimer water CO, Captive-Bubble 0.2-15 MPa; not significant changes water-wet 
CO,+N, ADSA 318 K with P 
Andrew et al. Carbonate Ketton KI, 7 wt% CO, Pore-Scale 10 MPa; 323 K weakly water-wet 
(2014) limestone (99% Calcite) microCT 
Al-Yaseri et al. Sandstone Quartz water, NaCl, CO, Nz Sessile-Drop 13 MPa; 333 K 47° 40.6° 33.9° 
(2015) 0.5 wt% CO)+N> 
Arif et al. (2017) Carbonate Calcite water, NaCl, CO, Tiltting-Plate 0.1-20 MPa; @ increases with P & strongly to weakly 
0-20 wt% 308-343 K salinity @ decreases water-wet 
with T & roughness 
Stevar et al. (2019) Carbonate Calcite NaHCO3 (aq), CO, Sessile-Drop up to 30 MPa; @ increases up to water-wet 
1 mol/kg Titling plate 298-373 K 10 MPa, after 10 MPa intermediate-wet 


@ decreases, and being 
constant @ decreases 
with T & roughness 


Haeri et al. (2020) Sandstone Navajo NaHCO3 (aq), CO, Captive-Bubble No clear trend with P & strongly to weakly 
Nugget Bentheimer 1 mol/kg T water-wet 
Bandera Brown Berea 
Mt. Simon 

Iglauer et al. (2021) Sandstone Quartz and NaCl, 10 wt% Hy Sessil-Drop 0.1-25 MPa, @ increases with P&T  water-wet less than 
aged Quartz, in acid Tilting plate 296-343 K 50°, for aged surface: 


intermediate-wet 


Table 10 Table 11 
Density of water and hydrogen phases used for Density of brine and hydrogen phases used for con- 
contact angle measurement of hydrogen/pure wa- tact angle measurement of hydrogen/brine (5000 ppm 
ter/Bentheimer tests. NaCl)/Bentheimer tests. 

Temp. (°C) Press. (bar) pw(kg/m?) Py, (kg/m?) Temp. (°C) Press. (bar) pą(kg/m?) Pu, (kg/m?) 

22.3 20.3 997.51 1.65 21.3 20.0 1000.91 1.63 

23.5 50.2 998.57 3.98 22.1 51.9 1002.16 4.13 

23.4 70.7 999.51 5.55 22.3 71.5 1002.99 5.63 

23.9 100.5 1000.71 7.73 22.9 100.5 1004.14 7.76 

24.4 23.5 997.16 1.89 24.0 20.3 1000.29 1.64 

24.5 50.7 998.35 4.01 24.8 49.4 1001.39 3.91 

25.0 70.4 999.11 5.49 24.7 70.9 1002.37 5.54 

25.1 100.7 1000.42 7.72 24.4 100.9 1003.77 7.75 

31.9 22.0 995.12 1.73 31.7 21.0 998.24 1.65 

32.5 51.8 996.25 3.99 32.2 49.9 999.36 3.85 

32.8 71.5 997.01 5.44 32.7 71.1 1000.13 5.41 

33.2 100.5 998.15 7.51 33.1 98.9 1001.21 7.40 

39.5 20.3 992.68 1.56 38.9 19.6 995.91 1.51 

39.9 50.2 993.83 3.78 39.5 50.8 997.05 3.83 

40.1 72.8 994.73 5.41 39.9 69.9 997.73 5.20 

40.3 100.3 995.84 7.33 40.1 100.1 998.95 7.32 

38.8 20.7 992.93 1.59 38.7 20.4 996.02 1.57 

39.0 48.8 994.08 3.69 39.0 51.1 997.24 3.86 

39.4 70.6 994.88 5.26 39.2 70.1 997.98 5.23 

39.5 99.2 996.07 7.27 39.4 100.4 999.20 7.36 

49.1 19.8 989.17 1.47 47.4 20.7 992.88 1.55 

49.2 50.6 990.44 3.70 48.3 51.3 993.82 3.76 

49.3 70.2 991.23 5.08 49.0 70.6 994.35 5.11 

49.3 101.2 992.53 7.20 49.2 100.7 995.53 7.16 
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Table 12 Table 16 
Density of brine and hydrogen phases used for contact Contact angle values of hydrogen/pure water/Bentheimer, repeated tests. 
angle measurement of hydrogen/brine (50,000 ppm 
NaCl)/Bentheimer tests. Test No. Temp. Press. Once Prange C) Vol. ave Vol. range 
(°c) (bar) (°) (mm?) (mm?) 
Temp. (°C) Press. (bar) pw(kg/m?) pn, (kg/m?) 1T-20°C 244 235 33.7 [30.5,38.8] 4.88 [1.99, 8.23] 
313 11 1004.18 166 245 50.7 348 [303,427] 4.05 [1.21, 7.10] 
319 514 100482 397 25 70.4 375  [33.4, 44.4] 2.75 [0.96, 4.41] 
330 706 100443 537 25.1 100.7 364 [32.4, 41.8] 3.91 [1.29, 8.37] 
333 100.7 100537 «7:52 2T~40°C 388 20.7 35.3 [315,428] 3.92 [1.30, 6.64] 
39 488 353 [32.4,41.6] 4.21 [1.41, 7.64] 
394 706 31.1 [27.3,349] 4.47 [1.79, 8.59] 
39.5 992 362 [31.6,42.3] 3.58 [1.15, 7.22] 
Table 13 
Density of brine and hydrogen phases used for con- Table 17 
tact angle measurement of hydrogen/synthetic seawa- Contact angle values of hydrogen/brine (5000 ppm NaCl)/Bentheimer, repeated 
ter/Bentheimer tests. tests. 
Temp. (°C) Press. (bar) —p,, (kg/m?) Pu, (kg/m?) Test No. Temp. Press. Oave Prange (°) Vol.ave Vol. range 
CO (bar) (°) (mm?) (mm?) 
31 20.5 1006.82 1.61 
31.3 50.2 1007.89 3.88 1 T~20°C 24 20.3 33.7 [29.7, 40.4] 3.92 [1.53, 5.56] 
31.6 69.6 1008.51 5.32 24.8 49.4 35.6 [30.8, 42.6] 3.63 [1.18, 6.79] 
31.6 100.9 1009.85 7.57 24.7 70.9 35.9 [30.5, 43.1] 3.45 [1.08, 5.97] 
244 1009 324 [295,375] 4.33 [1.54, 7.46] 
2T~40°C 38.7 20.4 31.3 [30.1, 33.3] 4.06 [2.41, 5.91] 
Table 14 39 51.1 31.69 27.4, 35.5] 3.89 [1.78, 6.27] 
Density of water and hydrogen phases used for con- 39.2 70.1 37.4 [34.5, 40.3] 2.54 [1.32, 3.90] 
394 1004 336 [305,386] 3.51 [1.18, 6.43] 
tact angle measurement of hydrogen/pure water/Berea 
tests. 
Table 18 
Temp. (°C) Press. (bar) — pw(kg/m?) py, (kg/m*) Contact angle values of hydrogen/pure water/Berea, repeated tests. 
23.6 20.8 997.23 1.68 
23.5 50.6 998.59 401 Test No. Temp. Press. O.,¢ Prange (°) Vol.ave  VOl-range 
23.7 70.2 999.42 5.50 Cg aban. 0) (mm?)  (mm°) 
23.9 100.7 1000.72 7.75 1T~20°C 24 20.7 30.5 [26.4,38.6] 5.17 [0.71,10.29] 
24.0 20.7 997.13 1.67 23.7 50.7 293  [25.3, 37.5] 5.35 (0.57, 12.05] 
23.7 50.7 998.55 4.02 235 71.5 266 [22.7, 32.8] 5.93 [0.76, 13.08] 
23.5 71.5 999.52 5.60 23.5 100.8 242 [22.2, 26.5] 9.68 [1.34, 17.42] 
23.5 100.8 1000.82 7.77 2T~20°C 224 19.7 228 [20.76, 25.06] 8.66 [2.16, 16.84] 
22.4 19.7 997.46 1.60 228 491 264 [23.44, 29.44] 6.02 [1.05, 11.97] 
22.8 49.1 998.69 3.91 
32.6 19.4 994.80 1.52 
32.7 50.0 996.11 3.85 Table 19 
32.8 69.3 996.92 5.28 Contact angle values of hydrogen/brine (5000 ppm NaCl)/Berea. 
33.0 101.1 998.24 7.56 
38.6 21.2 993.02 1.63 Test No. Temp. Press. Oavé Orange (°) Vol. ave Vol. range 
38.6 51.0 994.30 3.85 CO ba) © (mm?) (mm?) 
38.6 69.4 995.09 5.19 1 T~20°C 216 20.2 27.2 [23.3,30.9] 6.91 [2.86,11.48] 
38.9 100.7 996.33 7.39 218 506 28.7 ([25.3,31.8] 674 [1.41,15.52] 
47.6 20.5 989.78 1.53 223 70.2 296 ([25.6,33.3] 5.21 [0.97,12.11] 
478 49.4 990.93 3.63 226 1002 27.5 [253,301] 3.59 [2.15,5.03] 
48.2 70.6 991.68 5.12 2T~30°C 324 21.66 27.1 [251,283] 7.21 [3.05,11.91] 
48.2 99.7 992.90 7.12 324 506 27.8 [23.5,31.7] 5.99 [1.57,11.15] 
32.5 71.4 28 [23.7,31.2 4.19 [1.18,8.15] 
32.5 98.7 27.9 [26.0,30.9 5.05 [1.22,9.11] 
Table 15 4T~50°C 47.3 21.6 30.4 [26.5,33.2 4.71 [2.63,6.96] 
Density of brine and hydrogen phases used for con- 47.7 49.6 29.1 [26.6,33.7 4.3 [0.94,7.79] 
tact angle measurement of hydrogen/brine (5000 ppm 48.1 69.9 27.3 [19.9,34.5 3.88 [1.04,7.37] 
NaCl)/Berea tests. 485 996 29.3 [266,337] 4.86 [0.88,10.18] 
Temp. (°C) Press. (bar) pw(kg/m?) Pu, (kg/m?) 
21.6 20.2 1000.85 1.64 Declaration of Competing Interest 
21.8 50.6 1002.17 4.04 
aoe ee EAT ae The authors declare that they have no known competing financial 
32.4 21.6 998.06 1.69 interests or personal relationships that could have appeared to influence 
32.4 50.6 999.33 3.90 the work reported in this paper. 
32.5 71.4 1000.21 5.44 
ee eh es sp CRediT authorship contribution statement 
47.7 49.6 993.98 3.65 
48.1 69.9 994.68 5.07 Leila Hashemi: Methodology, Writing - original draft. Wuis 
48.5 99.6 995.77 7.10 Glerum: Formal analysis, Data curation. Rouhi Farajzadeh: Method- 


ology, Writing - review & editing, Supervision. Hadi Hajibeygi: Con- 
ceptualization, Methodology, Writing - review & editing, Supervision, 
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